On cosmologically designed modified gravity theories 



O 
(N 

o 

Q 



Jai-chan Hwang^, Hyerim Noh^, and Chan-Gyung Park^, 
^Department of Astronomy and Atmospheric Sciences, Kyungpook National University, Taegu, Korea 
^Korea Astronomy and Space Science Institute, Daejon, Korea 

Versions of parameterized pseudo-Newtonian gravity theories specially designed for cosmology 
have been introduced in recent cosmology literature. The modifications demand a zero-pressure 
fluid in the context of versions of modified Poisson-like equation with two different gravitational 
potentials. We consider such modifications in the context of relativistic gravity theories where the 
action is a general algebraic function of the scalar curvature, the scalar field, and the kinetic term 
of the field. In general it is not possible to isolate the zero-pressure fluid component simultaneously 
demanding a modification in the Poisson-like equation. Only in the small-scale limit we can realize 
some special forms of the attempted modifications. We address some loopholes in the possibility of 
showing non-Einstein gravity nature based on pseudo-Newtonian modifications in the cosmological 
context. We point out that future observations of gravitational weak lensing together with velocity 
perturbation can potentially test the validity of Einstein's gravity in cosmology context. 

PACS numbers: 98.80.-k, 98.80.Jk 
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I. INTRODUCTION 

In parallel with recent growth of cosmologically rele- 
vant observational data, theoretical cosmology also has 
been growing rapidly. Theories meet with observations 
mostly in the linearly perturbed stage of the large-scale 
cosmological structures. In our present understanding, 
the simple theoretical model of currently ACDM (cos- 
mological constant A as the dark energy together with 
the cold dark matter) dominated model, together with 
inflation generated initial fluctuations by a single field, 
has enough parameter space to amply cope with current 
observations. Recently, however, motivated by poten- 
tial future precision data, more complicated models for 
the dark energy sector have been pursued in the litera- 
ture. These include a fluid model with exotic equation 
of states, a scalar field with diverse potentials, modified 
relativistic gravity theories, etc; see Ref. [l| for recent 
reviews. 

One recent attempt concerns modifying perturbation 
equations by some unknown parameters in the level of 
modified gravity or pseudo-Newtonian gravity (2|-|4]| . Ef- 
fects of the parameters on the cosmic microwave back- 
ground radiation, the weak lensing shear field and the 
growth of large-scale structures, have been studied, and 
constraints on the parametrized modified gravity mod- 
els based on current and future observational data have 
been given. Proper theoretical justification of such at- 
tempts would be desirable. One way of justifying the 
situation is to show that such modifications are allowed 
or can be implemented in the relativistic gravity theories 
based on the action formulation. In this work we will ex- 
amine f{R,(j)^X) gravity as a potential candidate. The 
/(i?, 4>, X) gravity is presented as an action 



(1) 



where / is a general algebraic function of the scalar cur- 
vature R, a scalar field and the kinetic combination 
X = i(/)''^(/)_c; Lm is the matter Lagrangian. This in- 
cludes /(i?) gravity and many other gravity theories like 
Brans-Dicke theory, scalar-tensor theory, non-minimally 
coupled scalar field theories, etc, as cases. 



We will show that in general the pseudo-Newtonian 
modifications cannot be accommodated in our consid- 
ered generalized gravity theories, see Sec. [TTl Only in 
the small scales far inside the horizon we often have a 
closed form second-order differential equation of CDM 
(or baryon) density perturbation in a certain gauge con- 
dition with specific effective gravitational constant, see 
Sec. mil Such small scale coincidences in several mod- 
ified gravity theories may have motivated the pseudo- 
Newtonian approach. We point out some loopholes in 
the pseudo-Newtonian approach, see Sec. IIVI Even in 
Einstein's gravity theory with a minimally coupled scalar 
field as a dynamical dark energy, we can easily introduce 
the field potential where the dynamical nature of dark en- 
ergy is important in perturbation evolution. We believe 
it is desirable to explain the observational data based on 
theoretically motivated models, and also based on proper 
treatment of the complete set of equations provided in 
such models. In Sec. IIVI we also point out future obser- 
vations which can be used to test potentially Einstein 
gravity nature. 



The Appendix presents the complete set of equations 
in the f{R,(j),X) gravity in the presence of additional 
matter. These will be useful in future numerical study of 
perturbation evolution in such gravity theories. We set 
c = 1 EE 8nG. 
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II. VALIDITY OF COSMOLOGICALLY 
MODIFIED PSEUDO-NEWTONIAN GRAVITY 



Several versions of phenonienologically modified grav- 
ity theories specially designed for cosmological purpose 
have been introduced in the literature |2h51 iTHlll 

m, nn, m m mm m m, m, m [ii i^,iir The 

modified gravity assumes 



"fx 



where ip-^ — ip — Hx and 



(2) 
(3) 



a — X are gauge-invariant 
combinations which are the same as if and a, respec- 
tively, in the zero-shear gauge (x = 0); Sfiy = 5^ — fiav is 
a gauge- invariant combination which is the same as in 
the comoving gauge (v = 0). In Bardeen's notation [42| 
we have ip-^ = ^h, Q^x = and S^^ = fj,em- The gauge- 
invariant combinations and ip-^ can be interpreted as 
the Newtonian and the post-Newtonian potentials, re- 
spectively d^l- The variables 77 = f]{k,t) and jl = fi{k,t) 
are new free parameters introduced to modify Einstein's 
gravity where fj = 1 = p.. Later we will show that we 
could have p- ^ I even in Einstein's gravity. 

Our metric and energy-momentum tensor conventions 



ds^ = _ (1 + 2a) a^drf' - 2a^ 13 ^^dr^dx"' 
gS(l + 2^) + 27,„J dx'^dx^, 



(4) 



with X = 0/3 + 0^7, and 



A. Single component situation 

We consider a single component in an yet unspecified 
generalized gravity, and accept [1, [26l | 



(6) 
(7) 



where 5 = Sfi/fi. These equations follow from Eqs. (|50p 
and (pSj) , and Eq. , respectively, by demanding p ~ 
and 5p = = H. As we consider a single component, all 
fiuid quantities are about that single fluid supported by 
the nature of generalized gravity. In this case, because 
we assume 11 = 0, from Eq. (gl]), and Eqs. (gH]) and (|T7l) . 
respectively, we inevitably end up with 



-a 



X' 



'fx 

^ -Ax 
^2 ~ 2 



(8) 
(9) 



In the f(R,(j),X) gravity context we have /i = fix, 
p = px, Sfi = Sfix, Sp = Spx- {n + p)v = (/ix + Px)vx, 
and n = Tlx ] the subindex X indicates the contribu- 
tion from f{R,(j),X) gravity sector; the fluid quantities 
of the X-component are presented in Eqs. (|57p and (|62p . 
This apparently implies fj = 1 = p. That is, in the con- 
text of our generalized relativistic gravity theories the 
phenomenological modiflcations in Eqs. ^ and ([3]) with 
77 7^ 1 or /2 7^ 1 are not allowed. 

Even in the case such a fluid with j2 ~ 1 = fj is 
implemented in the context of f{R,(f>,X) gravity, thus 
/i = yUx a~~^ etc, those are at best uninteresting, 
because the nontrivial solutions in f{R,(f>,X) gravity 
should not affect the Einstein's gravity nature of the met- 
ric and fluid variables. 



= (p + 5p)5^p + 1 - i<5^nl^^) . (5) 

A vertical bar is a covariant derivative based on the co- 

(3) 

moving three-space metric g)^p of the Robertson- Walker 
spacetime. The complete sets of scalar-type perturba- 
tion equations in Einstein's gravity and in the f{R, </>, X) 
gravity are presented in the Appendix A and B, respec- 
tively. 

It is important to notice that, in the presence of mul- 
tiple component of fluids or a fluid in the context of gen- 
eralized relativistic gravity theories, the fluid quantities 
appearing in Eqs. ([3]) and ([5]) are collective ones, see Eqs. 
([53| in Einstein's gravity, and Eqs. ([6T|) and (|63)) in the 
f{R,4),X) gravity. In the following, except in the Ap- 
pendix A and B, we assume a flat background, thus set 
K = Q. 



B. Multi-component situation 

Now we consider a CDM component in an yet unspec- 
ifled generalized gravity theory, and accept 

k 

Scx = --Vcx-^'Px^ (10) 

Vcx + Hvcx = ^a^, (11) 

where a subindex c indicates the CDM component. 
Equations (HU]) and (HH) follow from Eqs. 1^ and P5|) . 
and Eq. (|55p . respectively, by demanding Pc = and 
(5pc = Etc. In this case, from Eq. (|48l), and Eqs. (|46| 
and (|47p . respectively, we have 

'^x = -«x-n, (12) 
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Thus, Eqs. ([2]) and ([3]) demand 11 = (77 — l)^^ and 



Later, in the small-scale limit, in Eq. (1271) wc will show 



1 

v' 



(14) 



This is apparently a stringent relation required between 
the two factors introduced in Eqs. ([2]) and ([3]). 

If a CDM component is present in the context of 
f{R,4>,X) gravity, from Eqs. (|5T|) and (p^ we have 
n = Ux = SF^/F; SF^ = 5F — Fx is a gauge-invariant 
combination which is the same as SF in the zero-shear 
guage (x = 0); we have F = df/dR. Thus, Eq. ([M]) 
demands 



SFx ^ in - 



(15) 



This is a strong physical condition on the nature of a 
gauge- invariant combination 6F^. 

More seriously, however, it is important to notice that 
Sfiv in Eq. ([T3|) is a collective density perturbation in the 
collective comoving gauge (v = 0); i.e., we have 



(16) 



where (5^ and (/i -\- p)v are sum of the fluid quanti- 
ties, see Eqs. (|53p and On the other hand, the 
CDM density perturbation in the CDM-comoving gauge 

is 5^cva = + 'i{a/k)H ^cVc = ^l^cx + 3(a/fc)7J/icfcx- 
If one of the components is due to generalized gravity 
theory like f{R,(t),X) gravity, the combination becomes 
more complicated. If we consider the CDM in f{R, cf>, X) 
gravity, from Eq. (pT|) we have 



{fi+p)v = — AicUc + ifJ-x +Px)vx, 
Sp^Spx, Ii = Ilx- 



Thus, we have 



r 



(17) 



(18) 



where 



5F,^=5F~F^vx. 



SfJ-x - fiXj^Vx, 



(19) 



Therefore, unless we have _F = 1 and Sfixvx ^ 
{jjLc/ F'^)5Fy^ = 0, both of which are quite nontrivial 
physical conditions on the nature of the generalized grav- 
ity, it is not possible to interpret 5^^ in Eq. (jSj as 5^,cvc 
which is 5^c in the CDM-comoving gauge (uc = 0). It 
is not likely that after implementing these conditions we 
have meaningful nontrivial results. For example, we can 
show that, in the f{R) gravity case, implementing -F = 1, 
SjjLx — if^c/ F^)^F = and vx = simply leads to an in- 
consistent equation for the background, thus impossible 
to implement. 



?7 - 1 1 



77 



1 F 



(20) 



Thus, compared with Eq. (fT8| . apparently the contribu- 
tion from X-component is important in our generalized 
gravity theory. 

In the context of generalized gravity theories, in gen- 
eral, it is unavoidable (except for the cosmological con- 
stant case) that we have additionally coupled second- 
order differential equation for the X-componcnt repre- 
senting the f{R,(p,X) gravity sector. The complete 
sets of equations of CDM perturbation in the context 
of f{R,(j),X) gravity, which inevitably lead to a fourth- 
order differential equation, are presented in the Appendix 
C. In the zero-shear gauge Eqs. ((TO)) and ([TT|) together 
with Eqs. ([70)1 and ((7T|) are the complete sets we need. 
In the Appendix C we also present the complete sets 
of equations in the uniform-field gauge and the CDM- 
comoving gauge. In all cases we end up having fourth- 
order differential equations. 



III. SMALL-SCALE LIMIT APPROXIMATION 



In the literature we often find a second-order differen- 
tial equation of Scv^ or ^cx derived in the small-scale limit 
[Tsl . [35I l36j . By naively preserving only terms attached 
with fc^ except for 6c term, Eq. (|55|) becomes 



2HSr- 



From Eqs. (^5)) and we have 



SF^ 



(21) 



(22) 



In the small-scale limit from Eq. ([T3| . using Eqs. ()17p . 
(I62ll and (HSl), we have 



fc2 
n2 



5Fy, 

'2F 



1 

2F 



(23) 



We need one more relation. From the equation of motion 
in Eq. §^ and using Eqs. and 5F = F^^Scj) + 

FjiSR [here wc consider F,x = case] we can derive 



SR^ 



1 

-2K 



f,chch + f. 



X- 



F.R A;2 



(24) 



where we have kept /^^ term compared with /.x^^/a^ 
term, and order 1 compared with {F^n/ F)k'^ / [see be- 
low Eq. (pS)) ] as in [l^lH, H^; we ignored f.x<pSX term 
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in an expansion of 5f,^ term in Eq. (|64p . From Eqs. 
I, we have 



1 + 2 — + J - 2—^^ 

^1 + 4— [J.<t><t> + J.x^) - 4— 
xa^ = -f]a^, (25) 
1 (^1 + 4— + -4— 

+ 3— [J.4><t> + Ix^) - ■5— 
xSficv,- (26) 



These can be compared with Eqs. ([2|) and However, 
it is important to notice that Sficv^ differs from Sfiy. By 
comparing Eqs. (IT3l) and (|23)) . and using Eqs. p2|) and 
(PSj). we have 



1 + 77 



(27) 



Using this, we have 
fc2 



1 1 



a2"^ - 1 + ?7F 



(5/ic 



Therefore, we have 



1 

V 



(28) 



(29) 



which is an aheady presented resuh in an exact context 
in Eq. ([T4)). The above relations are derived in the gen- 
eral context of f{R, 0, X) gravity which usually leads to 
higher order theory: see below Eq. ([56]). 
From Eqs. (HJ) and (UHl) we have 

r or.r 1 r 12/2 1^ 

Z Z 1 -j~ r 

= AttGcsShcv,, (30) 

where the effective gravitational constant becomes 

2/2 1 



SttG, 



cff 



(31) 



For / ~ R—2X — 2V we have a correct Einstein's gravity 
result with 87rGcff = 1. 

Now, we consider two cases: (i) / ~ RF{(j)) + 2p(0, X) 
theory, and (ii) / = f{R) theory. 

(i) In the case of / = RF{(j)) + 2p{(t>,X) theory, Eqs. 
(12^1) and (1261) reduce to 



^ + f4 



(32) 



2^ /,x-3i%)iU+/ 



-S^J.cv,■ (33) 



00 



(ii) In the case of / = /(-R) theory, however, the equa- 
tion of motion in Eq. (j64p is not available. Instead, from 
Eq. (|65|) we have 



*a2 



1 



^^x - ("x + 2<^x) - f^tJ-cx " „2 



a2 F 



(34) 



We take Sex 
(I34l) and (IM 



JctJc in the small-scale limit. From Eqs. 
we recover Eq. ([23)) . Now, in the f{R) 
theory we have SF = F fjSR. Using this relation, from 
Eqs. (I22|) and ([23]) together with Eq. ([Ml) we can derive 



"Px 



1 1 



2^^1 + 3^ 



5^c 



(35) 
(36) 



Although these results can be regarded as cases of Eqs. 
([25l) and ([26l) . we need a separate treatment because Eqs. 

and (|26p were derived from the equation of motion 
which is not available in f{R) gravity. As in Eqs. ([M)) 
and (pS]) we kept (5i? = SF/F^n term even compared with 
[k"^ / a?')5F / F term because becomes zero in the Ein- 
stein gravity limit. Thus, it is important to notice that 
even in the small-scale limit we are already considering, 
we have to set F^^/F <^ a? / <Si 1 to properly have the 
Einstein's gravity limit. 

Apparently, the derivation of Eq. ([50)) in the small- 
scale limit requires rather complicated steps, and the only 
way to check its validity is comparing with the exact re- 
sult. In Fig. [T] we estimate the accuracy of the above 
asymptotic approximations using a case of f{R) gravity. 
We consider a case with / = R^^" + qR-'\ The R^^" 
term dominates in the early (radiation) era and allows 
a scaling evolution of and Sfix following the domi- 
nant fluid; first follows radiation and then matter. The 
qR~'^ term can be tuned to cause accelerated evolution 
in recent era; we have studied this case in |44[ | in de- 
tails. For the perturbations we solve equations in the 
CDM-comoving gauge presented in Eqs. ([551) - ([55)) with 
additional presence of radiation component. 

The Figure shows that in the small scale far inside 
the horizon the approximate equations in Eqs. (j2ip and 
([36)) reproduce the evolution of Scv^ properly. That is, 
as the scale becomes smaller results from the small-scale 
approximation more properly reproduce ones from full 
calculation. The evolutions of Gcff, however, show that 
deviation of Gcfr from Newtonian G is more significant 
as the scale becomes smaller, see Eq. ([55)) . In practical 
applications it is always important to check the validity 
of such asymptotic form approximations compared with 
exact solutions based on the complete set of equations. 
The complete set of perturbation equations in several 
gauges in our generalized gravity theories is presented in 
the Appendix C. 
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FIG. 1: Top panels from left to right: Evolution of density fn {i = r,m,X); recent behavior of density parameters due to 
f{R) gravity; test under the type-la supernovae observation; evolution of the Hubble parameter. We use a functional form 
f{R) = R^'^'^ +qR~^ with e = 10~^ and n = (red), —0.1 (yellow), —0.2 (blue curves). Note that we compare /ir and fim with 
Fnx- Evolution of F^x shows early scaling and its recent domination as a dark energy. We use definitions, Q.i = fij / {3F H^) 
for radiation and matter and Q,x = fJ-x /{SH'^). The binned type-la supernovae data are based on the Union2 sample [4^, and 
Hubble parameteter data are from Ref. [i^. Bottom panels from left to right: Behaviors of the CDM density perturbation 
growth factor in three different scales; the last figure shows evolutions of Gcff. In the growth factor panels, results from the 
full calculation (including radiation; solid curves) and from the small scale approximation (dashed curves) are shown together. 
Green dots with 1% error bars are expectations from the future X-ray and weak lensing observations As the scale becomes 
smaller results from the small-scale approximation more properly reproduce ones from full calculation. The evolutions of GcS, 
however, show that deviation of Gcff from Newtonian G is more significant as the scale becomes smaller. 



IV. POTENTIAL COSMOLOGICAL 
SIGNATURES OF NON-EINSTEIN GRAVITY 

Can we find evidence of non-Einstein gravity nature 
from the large scale cosmological observations? Be- 
fore presenting a (in principle) possible way of using 
cosmological observations to distinguish potential non- 
Einstein' gravity nature, here, we would like to point out 
some potential loopholes in recent such arguments based 
on pscudo-Ncwtonian approaches motivated by approxi- 
mate treatment of modified gravity theories. 

One argument is based on the deviation of baryon (or 
CDM) density perturbation equation and Poisson's equa- 
tion in modified gravity. In the small-scale limit from 
Eqs. dSI]), (1201) and §^ we have 

6cv, + "^HScv, — — ^Qfx, (37) 
-^^x ^ -inGcffS^icv,- (38) 

Under approximations used in previous section we have 
Eq. ()31|) . and in the Einstein's gravity limit we have 
87rGoff = 1. Here, it is important to remember that 
the closed form second-order equation in the considered 
modified gravity theories is valid only in the small-scale 
limit. 

Although in Fig. [1] we have shown an example where 
the CDM (or baryon) growth rate is well approximated 
by the small-scale limit approximation in Eq. (j30|) . that 



is not necessarily the case for other cosmological observa- 
tions like the CMB and the density power spectra, etc. It 
is important to notice that even in Einstein's gravity, in 
the presence of dynamic dark energy, in general, we can- 
not ignore the contribution of the dark energy perturba- 
tion in the right-hand-side of the Poisson's equation: see 
Eqs. pT)) and One such an example demonstrating 

the importance of coupling even in Einstein's gravity is 
presented in (isj . 

In [i^ we have shown that even in the case of Ein- 
stein's gravity we have observationally significant devi- 
ations by ignoring the dark energy perturbations in the 
CMB and density power spectra, and in the baryon den- 
sity perturbation growth rate. As we emphasized in this 
work, consideration of dynamical dark energy in general 
causes a coupling between the baryon (or CDM) den- 
sity perturbation and the dark energy perturbation; and 
the dynamics of dark energy perturbation is described by 
its own additional equation of motion which is again cou- 
pled through gravity with other components. This is true 
even in the context of Einstein's gravity with a minimally 
coupled scalar field as a dynamical dark energy. 

In [4^ we have studied a case of dark energy model 
based on a minimally coupled scalar field with a double 
exponential potential as an example. We showed that 
proper consideration of the dark energy perturbations is 
important in comparing theories with present day cosmo- 
logical observations. In the CMB temperature anisotropy 
power spectrum, the large-scale density power spectrum, 
and even in the perturbation growth rate we have shown 
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FIG. 2: Top panels: Evolution of SnGcS as a function of 
scale factor a{t) and the redshift z for several different scales. 
As the scale becomes smaller SnGeS approaches unity. We 
consider a double exponential potential V{(fi) ~ Fie'^^** + 
V2e~^^'* where (j) is the scalar field and A2 = 1.0; for other 
parameters, see [3. Bottom panels: Evolution of the ratio 
of the scalar field compared with the ordinary fluid (matter 
plus radiation) using the same model as in the top panel. 
Deviations of the ratio from unity in the early radiation and 
matter dominated eras are due to the scaling nature of scalar 
field we considered. 



that it is critically important to take into account of 
dark energy perturbations properly. This result can be 
regarded as an evidence that the coupling between the 
baryon (or CDM) density perturbation and the dynami- 
cal dark energy perturbation is in general important even 
in Einstein's gravity. In the presence of a scalar field dark 
energy, in the CCG, from Eqs. ([M ]) . ([55 ]) and (P ]) 

we have 



6c + 2HSc = 



F 36p) 
A-KGeeSfJ-c- 



- 4(1)6(1) - 2V^S<p 



(39) 



In Fig. [2] we present evolution of STrGoff and the effect of 
the scalar field contribution using a double exponential 
potential model of a minimally coupled scalar field in 
Einstein's gravity |48|. Although becomes negligible as 



we go to the small scale, we have STrGoff 7^ 1 even for 
a minimally coupled scalar field as the dark energy in 
Einstein's gravity. 

The other argument is based on Eq. so that fj ^ 
1 could be regarded as a strong signature of the non- 
Einstein gravity nature. The main point used in this 



argument is that the gravitational weak lensing effect can 
measure the sum of Newtonian potential (— ct^) ^nd the 
post-Newtonian potential {}Px)- From Eq. (|48l) we have 



+ = -n, 



(40) 



Although we could have significant amount of anisotropic 
stress n even in Einstein's gravity, that could be indeed 
regarded as highly peculiar situation almost comparable 
to rather considering the non-Einstein gravity. As we 
have ~ —a^ in Einstein's gravity, we can accept that 
it is reasonable to regard 11 = {fj—l)a-^ ^ 0, thus 7? 7^ 1 as 
a strong signature of non-Einstein's gravity. Now, as a^ — 
can be measured from the weak lensing observation, 
the remaining problem is whether we can determine 
or ipy^ individually in Einstein's gravity. 

Here, it is interesting to point out that the evolution 
of baryonic velocity perturbation v^xi in principle, can 



determine the Newtonian potential 



|34| . Indepen- 



dently of the gravity theories, from Eq. (|55|) or (pT]) for 
the baryon component, we have 



(41) 



Similar equation for the CDM component is presented 
in Eq. (jlip . Thus, apparently, the velocity perturbation 
depends only on the Newtonian gravitational potential 
(Xx- Since this equation is derived from the momentum 
conservation of the baryon component, unless there exist 
exotic interaction with other components, this equation is 
generally valid independently of the nature of the gravity 
theory based on Riemannian geometry. 

Therefore, in Einstein's gravity theory we have 



(42) 



In this expression the left-hand-side and the right-hand- 
side can be determined from weak lensing observation 
and peculiar velocity field observation, respectively. By 
combining these two informations we can, in principle, 
tell whether the post-Newtonian gravitational potential 
(}Px) has the same amplitude as the Newtonian poten- 
tial i^—OLx)- If the two amplitudes differ from each other, 
i.e. if Eq. is violated, assuming validity of linear 
theory, that can be regarded as an important signature 
of non-Einstein's gravity nature revealed in the cosmo- 
logical situation. For recent study in this direction, see 



V. DISCUSSION 

We have considered feasibility of modifying Newton's 
gravity in cosmology from the perspective of generalized 
gravity theories. As a concrete example we examined 
/(i?, 0, X) gravity which includes /(i?) gravity and many 
other gravity theories known in the literature as cases. 
For the baryon or CDM component in the context of the 
generalized gravity theory, the perturbation equations 
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are inevitably coupled with the presence of the general- 
ized gravity theory sector, thus in general ending up with 
at least fourth-order differential equation; situation is the 
same even in the minimally coupled scalar field as a dark 
energy in Einstein's gravity. The closed form of CDM 
density perturbation equation proposed in the literature 
is only possible in the small-scale limit in the context 
of generalized rclativistic gravity theories. Even in such 
a small-scale limit we can easily introduce cases where 
proper treatment of generalized gravity is important to 
compare with the CMB and density power spectra. 

It is always important to analyze properly (i.e., with- 
out approximation) the roles of dynamic dark energy 
perturbation in Einstein's gravity before one concludes 
against the canonical gravity theory which has been ob- 
scrvationally successful in all weak gravity tests. We 
also have addressed potential loopholes in arguing non- 
Einstein gravity nature based on cosmologically modified 
pseudo-Newtonian gravity theories. We point out that 
future precise observations of the weak lensing together 
with the baryon velocity perturbation can potentially test 
non-Einstein gravity nature based on cosmological obser- 
vations. 



where x = cip + ct^7 is a spatially gauge-invariant combi- 
nation. We note that all our perturbation equations are 
spatially gauge invariant. Although derived in Einstein's 
gravity these equations are valid even in the context 
of generalized gravity theories including the f{R,(t),X) 
gravity presented in the Appendix B. The generalized 
gravity theory we consider can be written in the Ein- 
stein's gravity form 



Gab -— Ta, 



(52) 



Appendix 

Appendix A: Perturbations in relativistic gravity 

In the Appendix A and B we consider the presence of 
general K. Background evolution is described by 



1 K 

If^ 2' 



/i — —3H + p) . 



(43) 



3 a' 

In the multiple component case, we have fi = and 
p = pj , together with 

-3Hi^M+p,) + Io^, (44) 



where Iqj 0. Based on our metric and energy- 
momentum tensor conventions in Eqs. ^ and (O, the 
scalar-type perturbations are described by [49j 



where Tab is reinterpreted as the effective energy- 
momentum tensor [50j . In this way, as long as the gener- 
alized gravity can be expressed as in Eq. (|52l) , the above 
basic set of perturbation equations is valid in those grav- 
ity theories as well. 



In the presence of multiple component fluids (including 
dust, dark matter, radiation, fields, etc.), we have 



K — —3(p + 3Ha H — -X, 



1 



-|- Hk = 



fc2 - 3K 



fc2 - 3A' 
K 5 — X 



3 , .a 
-(/i+p)-w, 



X + - - a 



8ii + 3H ((5/i + 5p) = (/i + p) K - 3Ha v 

' a 



(45) 

(46) 

(47) 
(48) 

(49) 



(50) 



\a^{[i^p)v\ 
a^ili+p) 



k 

-a + 
a 



a{p. +p) \ 3 



(51) 



i j 

{^l + p)v = ^(mj + Pj)vj , n = ^ Bj . (53) 



In the case of generalized gravity theory, see Eqs. (|5T|) and 
(j63p . The individual component follows its own equation 
of motion 
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+ SH {Sfi^ + 5pt) {fit +pt) ( K - 3Ha Vi] Woi, (54) 

\ a J a 

[aH,.+pM]- _ ^ k_ / ^ _ 2 A^^^ _ ^ ^ ^^^^ 



where J2i ^hj = = X), 51^. 

r 



Appendix B: Perturbations in f{R,(f>,X) gravity 



effective ones 



We consider /(i?, </>, X) gravity considered in |50|. The 
action is 



5* 



f{R,(f>,X) + L„ 



(56) 



where X = The minimaUy coupled scalar field is 

a case with f = R — 2X — 2V{(p). Even in the absence of 
additional matter component the f{R,(j),X) gravity, in 
general, leads to a fourth-order differential equation for 
the scalar-type perturbations. In order to have second- 
order differential equation with nontrivial F = df/dR 
we should have either (i) / = RF{(t),X) + 2p{(j),X) or 
(ii) / = /(i?),see 0. 

The gravitational field equation following from the 
above action can be arranged in the form of Einstein's 
gravity in which the new contributions are interpreted 
as effective energy-momentum tensor; the effective fluid 
quantities are introduced in Eq. ([5])- Using this strategy, 
the background and perturbation equations in Einstein's 
gravity remain valid with the fluid quantities replaced by 
the effective ones [sH - 

Background evolution is still described by Eq. 
with the fluid quantities reinterpreted as the following 



M = 

MX 



1 



^■X, P=—Pr, 



PX, 



PX^J 



1 / FR-f 



F + 2HF 



(57) 



where X = —hi^- The equations of motion is presented 
in Eq. (80) of [5| 



1 



0. 



We also have 



R = 6[2H 



■ K\ 

H+—]=fi-3p. 



(58) 



(59) 



The ordinary matter part follows its own equation of mo- 
tion /i„i = —3H (/i„i + Pm)- Thus, we have 



F 

fix + iH {fix + Px) = -^^^r. 



(60) 



J 



In our generalized gravity context Eqs. ()45p -(|5i p remain valid with the fluid quantities reinterpreted as the following 
effective ones 



1 , , 6F ^ 1 ^ 

Ofi = -p^Oflm + OfLx - fhn-p^, Op = —dp„i + Opx - P, 
{fL +p) v = ^ {flm + Pm) V,n + (MX + Px) VX, II = -^II,, 



n 



X, 



(61) 



with 



- {IxSX - /, 



Sp 



X 



1 

F 



1 



k (\ 



aF \ 2 



SF 

+ df^xX ~ 3HSF +[^f- f,xX + 3HF -F'—]—+Fk + 3HFa 



{f,xSX + /,^(50) + 5F + 2H6F - [ 7:f + F + 2HF - F 



1 



2fc2\ SF 2 • 



3a2 



^ -^FK-Fa-2(^F + HFy 



.X 



-^5F + HSF + Fa 



n 



X 



F 



5F- Fx] , 



(62) 
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and 5X 



+ (fP'a. If wc have multiple fluid components, we have 

3 3 

Sp^rn=^6fij, 6p.,n=^5pj, {pm. + P7n) V,n ^ ^ {fJ-j + Pj) Vj , Ilm='^llj, 



(63) 



where the individual component satisfies Eq. (j44p for the background and Eqs. (IMl) and ([55]) for the perturbations. 
The equation of motion is presented in Eq. (81) of [5^ 



f,x 

We also have 



/. 



X 



+ 2/,^a + 4 fa' ) +<5/,0 = O. (64) 



-/i - 4i7K + ( — - 3ff I a + 2 ^ 



(5i? = 2 

If the ordinary matter part follows its own equation of motion 



= 5jjL — 3Sp. 



Slim + 3H {Sfirn + Spm) = {pi,n + Pm) \ K - ^Ha - -V„i ] , 



rn 

~ —a 



Sp„ 



a'^il^rn+Pm) a a{flrn + Prn) \ 3 a'^ 

which follow from Eqs. ([54| and (f55|) . from Eqs. ([50]) and ([5T|) we have 



-n„ 



(65) 



(66) 
(67) 



6fj,x + 3H {6fj.x + 6px) - (px +Px) ( - 3Ha - ^vxj = -^Vm + ( SF - '2j;Sf \ , (68) 



^ [a^(/ix +px)vx] - - 



ifJ-x +Px) a + Spx - - 



2 fc2 - 3K 



3 a2 



-n 



X 



F ^ 

= iPm + Pm) Vm - ^^Pni'^^- (69) 



1 k 

F^a^ 



Notice that in the single component case without m-component. Eqs. (|55)) and (|69p are the same as Eqs. and 
(p7)) with the sub- indices m replaced by X. 



Appendix C: CDM perturbation in f{R,(fi,X) gravity 

Here we set K = Q. For simplicity, we consider (i) / = RF{4)) + 2p(0, X) and (ii) / = f{R). In the following, 
we present complete sets of equations of CDM in the above generalized gravity theories, in three different gauge 
conditions. These are the zero-shear gauge, the uniform-field (or F) gauge, and the CDM-comoving gauge. 



Appendix C.l: Zero-shear gauge 

The zero-shear gauge takes x = as the temporal gauge condition. Equations (|10p and (|TT|) are valid for differential 
equations for 5cx and Vcx- Using Eq. (jl2p can be expressed in terms of (p^ and SF^- Thus, two additional first-order 
differential equations for Lp^ and SF^ will complete the perturbation equations. One equation follows from Eqs. (j45p 
and (HTI 



,x 



, • 1 / , F\5F a 

SF \ H -\ Mc-Wc 

2 I F F ^"k 



(70) 



10 



The other one can be derived from Eq. p3 



1 • • 3F ■ 1 

^{f,x + 2f,xxX)cl,Scl,~—SF+- 



n2 



SF 



'2F 



(71) 



where we have SF = F,^54) for (i), and 5(t) ~ Q = X for (ii). From these we have the two equatfons for ip-^ and 5F^. 
And these two equations are the additional equations we need to solve together with Eqs. ((TU)) and PT|) . 



Appendix C.2: Uniform-field gauge 



We take the uniform- field gauge for (i), and the 
uniform- gauge for (ii); in both cases we set 5(j} ~ 
= 5F as the temporal gauge condition, and we call 
the gauge condition as UFG. We have 5X = —2Xa and 
5f x = -2f,xxXa. From Eqs. (gSJ-gZ]), Eqs. (gS]), gT]), 
and (|48p . respectively, we have 



F + 2HF 



^ + f,xX + 2f^xxX^ 

frr F \ 3a 



— 



2F 



5fic 



1 a 
2Fk 



fJ-cVc, 



aF 



aF 



+ 2F 



^ + f.xX 
F + 2HF 



Mc 1 / a 



(72) 



(73) 



where from Eqs. (P5 |) - (|T7)) we have 
1 



i la I 
2Fk^ ' ' 

Slic 



-±,(3HF- f.xX-2f^xxX' 



3 a 



3i^ 



(77) 



(78) 
(79) 



These equations give a set of differential equations for $, 
(5c and Vc] these can be combined to give a fourth- 
order differential equation. $ and 4* are gauge-invariant 
combinations, and all the other perturbation variables 
are evaluated in the UFG; 5c and Vc evaluated in the UFG 
are equivalent to gauge-invariant combinations S^sf = 
5c + {3HF)5F and VcSF = i'c — i^/aF)5F, respectively. 
The CDM density perturbation in the CDM-comoving 
gauge can be constructed by 



where 



5cv^ = 5c + SH-Vc, (80) 
where the right-hand-side can be evaluated in the UFG. 



In the context of /(i?) gravity we simply set X = 0. In 
the absence of the GDM component, these equations are 
presented in Eqs. (85) and (86) of [H. For the CDM 
part, from Eq. ([54|) and Eq. (|55|) . respectively, we have 



5c = Vc + K — 3Ha, 

a 

Vc + tlvc = -a, 
a 



(75) 
(76) 



Appendix C.3: CDM-comoving gauge 

The CDM-comoving gauge (CCG) takes = as the 
temporal gauge condition. From Eq. (|55p we have a = 0. 
Equations jSH) and (gT]), Eq. (08]), and Eqs. (|45]) and Eq. 
P7)) . respectively, give 



5c 



(h. 



1 

2i^ V 2 



-J.x 



F 
F 



X = ^ + 



6F 



-5F + H6Fj . 



5F + H6Fj , (81) 
(82) 
(83) 
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In order to complete we need a first-order differential 
equation for 5F. From Eqs. P5|) and (|Tf)) we have 



F 



3F • 1 

dF+ - 

2F 2 



- 2/,^xX + 3 if + — 



H 



1 
F 



-J - f.xX 



fc2 
^2 



fc2 

,5F= -2i^:3</j + (5/i,+ (i^ + 2i?i^) -2%, (84) 



where we have 5F = F,cj,S4> for (i), and S4> = = X 
for (ii). Equations ([5T|) -(|M | provide a complete set of 
equations for Sc, Xi V ^'Hd 5F. 

Equations ([5T|) - ([55)) can be combined to give 



5c^ is 5c in the uniform-curvature gauge {ip = 0) and ipg^ 
is the spatial curvature perturbation tp in the uniform- 
CDM-density gauge {6c = 0). 



1 



{6c + Sep) 



(85) 



where -I- 3(/3 = (5cip = 3ips^ is a gauge-invariant com- 
bination, thus can be evaluated in any gauge condition; 



In the case of f{R) gravity we suggest the following 
alternative form. From SR = Sfi— 3Sp, Eq. and Eq. 

(1511) we have 



SF + 3H5F +(-Ir+^)5F+ \f5R = Fk+\ {S^^n - 5Spm) , (86) 
\ 3 a^/ 3 3 

/ f\ SF ■ P\ SF ^ 1 

k+i2H--\K = -3H— + ^-f + mF - - F- j —--5R+ -5^im, (87) 

he K. (88) 
For a CDM component we have = <^Mni = Mc<5cj and ( 5pni — 0- 
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